Interferometric Parallax: 
A Method for Measurement of Astronomical Distances 
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We show that distances of objects at cosmological distances can be measured directly using in- 
terferometry. Our approach to interferometric parallax comes from analysis of 4-point amplitude 
and intensity correlations that can be generated from pairs of well-separated detectors. The base- 
line required to measure cosmological distances of Gigaparsec order are within the reach of the 
next generation of space-borne detectors. The semi-classical interpretation of intensity correlations 
uses a notion of a single photon taking two paths simultaneously. Semi-classically a single photon 
can simultaneously enter four detectors separated by an astronomical unit, developing correlations 
feasible to measure with current technology. 
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There is no more important problem in astronomy than 
resolving the third dimension of source distances. The 
crisis of dark energy and dark matter in cosmology hinges 
on distance measurements. Estimates using red shifts 
or Type la supernova sources are weakened by model 
dependence of the cosmology and assumptions on the 
evolution of distant sources. Here we show that direct 
measurements of cosmologically distance objects can be 
made from the analysis of correlations of detectors sep- 
arated on the scale of the solar system. Correlations of 
amplitudes (first order coherence) are used in Michelson's 
interferometer and radio telescopes. The breakthrough 
of Hanbury-Brown and Twiss demonstrated 
2-point correlations of intensity (second order coherence) 
developed by counting photon fluxes in separated detec- 
tors. Higher order correlations have been proposed ear- 
lier for reconstructing the phase of the coherence function 
2] and to improve the sensitivity 3] in intensity corre- 
lations. We show that 4-point amplitude and intensity 
correlations contain further information on the distance 
to the sources. The baseline required to measure cosmo- 
logical distances of Gigaparsec order are within the reach 
of the next generation of space-borne detectors. Measur- 
ing source distances of Megaparsec order appears feasible 
now. 

Let Xi be the position vectors of the ith detector rel- 
ative to the origin (Fig. QJ. Consider a point source at 
r — fr relative to the origin, also located at position rj 
relative to each detector: 

fj = r- xj. 

Calculate the distance from the Jth detector to the 
source to order 1/r: 



r.j = r — f ■ xj + 
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FIG. 1: Four detectors are positioned so that the relative 
position vector Aa?i2 between one pair is approximately the 
same as between another pair Aa;34, with a net translation X 
between the pairs. 



r and will be responsible for probing it via interferomet- 
ric parallax. The frequency-domain Green function for 
propagation from source to receiver is expanded 
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where 8^ (r) = 5 iJ — f l f J . Here upper indices denote vec- 
tor components. The third term depends on the distance 



where k is the wave number. We drop the e tkr /r pref- 
actor which cancels out in calculations. The formula ex- 
tends trivially to two distinct sources, S and S' , for which 
primed symbols such as f = f'r' take the obvious mean- 
ing. 

Concentrate for a moment on two detectors 1,2. Stan- 
dard analysis of the interferometric correlation treats the 
electric field E as a random variable, described by con- 
volution of the Green function with correlations of the 
source. Then the amplitude correlation for a single po- 
larization between the two receivers can be expressed as, 
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where Eq is the field of one source, Iq =< EqEq >, and 
so on for primed variables. Here 
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The overall phase ip — f ■ [xx — X2) — {x\x\ — x 2 x 3 2 )S l j!, /2r. 
The overall phase cancels in a determination of the abso- 
lute value of the correlation. In eq. the source-detector 
distances are consistently replaced by r and r' everywhere 
except in the phases. 

The intensity-intensity correlations, < Ixl 2 >=< 
\E(xi)\ 2 \E(x2)\ 2 > between receivers 1 and 2 for clas- 
sical light are found to be, 
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Here 3? denotes the real part. A quantum mechanical cal- 
culation gives a similar result, and incorporates photon 
bunching. We see that ( j ) P araUax appears in both am- 
plitude and intensity correlations. Intensity correlations 
can be used at optical frequencies by simply counting 
photons, and have certain advantages in automatically 
canceling the overall phase ip. 

The first phase (j)^ has been used extensively for an- 
gular position measurements. Notice that this phase is 
translationally invariant - it depends only on the relative 
position vector of the two detectors. The scale Axi2 is 
conjugate to a difference of wave numbers, not the ra- 
diation wavelength. Our focus is on the second phase 
^parallax ^ c onsu iti n g Fig. ^ let detectors 1 and 2 be 
located at average position X and separated by A212: 
x\ = X — Ai?i2/2, X2 — X + AX12/2. By substitution 
the second phase is 
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This phase will change with X even if Aa?i2 is fixed. The 
explanation of course is parallax sensed via curved wave 
fronts. A measurement of the correlation's dependence 
on translating a fixed detector pair in X can probe the 
distance dependence on 1/r and l/r'. 

We emphasize that interferometric parallax is quali- 
tatively different from standard trigonometric parallax 



long used in astronomy. In trigonometric parallax the 
distance is estimated by a precise measurement of the 
angular position of the source from two different loca- 
tions. For interferometric parallax a precise measurement 
of the angular position is not required, and dependence 
on the angular position is weak. Strictly speaking, the 
individual sources need not be resolved. Sensitivity exists 
in selecting source for the measurement, and excluding 
background. 

We now discuss finite source size effects. The basic 
two point amplitude correlation for an incoherent source 
component located at r. \ + y is 

< EtfjETfa) > = Jd 2 y< E Q (y)E*(y) > 
p ik[\r'i+^\-\r2+i/\ 
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Here we have expanded the argument of the exponential 
integrand as 

\n +y\- |r 2 + y\ = n - r 2 + fx ■ y - f 2 ■ y 
-(y-rx) 2 /2r+(y-f 2 ) 2 /2r. 

By inspection the last two terms only contribute at or- 
der l/r 2 . Let a be the effective transverse size of the 
source. All terms involved in X are functions of the an- 
gular size of the source a/r, and X is negligibly small if 
a/r >> l/(kAx). As long as Ax < a/(kr) the correc- 
tions due to finite size can be absorbed into the over- 
all factors Jo, Iq etc. in Eq. and a source appears 
to be a point. This reproduces the usual planar source 
criterion 4| that an observable signal requires baselines 
smaller than the coherence zone of each source. All 
higher order correlations for two sources can be expressed 
as products of such two point amplitude correlations of 
individual sources, justifying the point source approxi- 
mation. 

Orders of Magnitude: For wavelength A and typ- 
ical source angular separation A9 we have (f)^ ~ 
A0/(A/Axi 2 ) and ^ aUax ~ X/r/(X/Ax 12 ). We rec- 
ognize A/AX12 as the lower limit on angular resolution 
from optics. Similarly, the parallax phase is of order one 
if the baseline of translation X could be resolved by an in- 
strument of aperture Ax looking from distance r. Using 
the lower limit of the single-source coherence kAx < r/a 

iparallax fvt \/i / /\ iparallax j 

gives (p\ 2 ~ (X/a)(l — r/r ). For (p'{ 2 ~ 1 and 
comparable r, r' the source-size should match the trans- 
lational scale. Of course the coherence zone criteria do 
not require literally small sources, but represent the ex- 
istence of Fourier modes (structure) in the regime of size 
indicated. Setting <$^ a ax ~ 1 yields the distance scale 
that can be observed: 



r < 1 Gpc 
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Although phases can often be measured with exquisite 
accuracy, we will continue assuming (j)P arallax ^ 1 f or our 
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estimates. Consider detectors separated by 10 4 km, a 
near Earth orbit, and translating over X ~ 1AU in a 
period of a year. Numerically 
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A# 1mm Axi2 
arcsec A 10 4 fcm' 



.parallax , n -4 ^ lGpC 1mm Axi 2 
^ 2 ~ 10 
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Baselines of order 10 4 km at cm wavelengths have been 
demonstrated by current technology. For Gpc dis- 
tances, one needs to measure a relatively small phase 
or push the limits of baselines to much larger than 10 4 
Km and/or wavelength to the sub- mm range. This 
may be possible due to the huge range of possibilities 
for AO. Quasar sources are believed to have phys- 
ical sizes extending to the range of 1 AU, whereby 
AO ~ 10~ 9 arcsec at distances of Gpc order. The max- 
imum coherence zone for such sources Ax ~ Xr/a ~ 
10 n m (\/mm)(r/Gpc)(AU/a) are compatible with base- 
lines of order AU. Black hole and GRB sources are of 
course even smaller, with correspondingly larger coher- 
ence zones. It might also be possible to measure gravita- 
tionally lensed single objects, exploiting two path lengths 

/ t • • i j_ r i parallax r j 

r, r . In principle measurement ot tp\ 2 ot order unity 
can measure distances to Gpc order provided suitable 
sources can be exploited. 

The ratio <f><® j^raiiax _ 1Q 8 for a ty p ica i va i ue f 

AO ~ O.larcsec assuming Gpc distance. There are rea- 
sons to expect high control over by technological 
means. However for the rest of the paper we consider 
"worst case" scenarios, in which control of </>(°) is less 
than ideal. There happens to be a practical strategy 
to null out the effects of the rapidly varying phase. To 
control the effects of 4>^ we can (in effect) measure it 
twice, using another pair of detectors 3, 4, separated by 
the "same" offset: AX34 = Ax*i2 + if ■ It is clear that 
ff « Axi2 can be made relatively small with great pre- 
cision. We consider the product of < hh >< >, 
which is one of the terms in the 4-point intensity correla- 
tion. Denoting the overall normalizations by Ni and N2 
we write 

< hh >= Ni + N 2 cos(A + B); 
< hh >= N{ + N' 2 cos(A' + B'), 



where A 



B — <fi 1 y 2 iraUax , with the primes switch- 
ing 12 — > 34. In an experiment where A >> B and 
Ax\2 ■ (r — f') varies rapidly the products with an odd 
number of cosines will average to zero. The average here 
(symbol << >>) might occur over running time in which 
drifts of the detector position values cause <f>^ to vary. 
Another cause of variation lies in small ranges in u> differ- 
ing between the detectors, and there are no doubt other 
possible causes. The term with two cosines gives 

<< cos(A + B)cos(A' + B') »= 
< < cos A cos A' » cos B cos B' 
+ « sin A sin A' >> sinBsinB'. 



The only non-zero data will come from << 
cosAcosA' >>=<< sinylsinA' >>= 1/2, namely 
those regimes when the rapidly varying terms coincide. 
Collecting the terms gives 

< hh >< hh >= N X N{ + N 2 N 2 cos(B - B') 

= NiN[ + N 2 N 2 cos(^ rallax - cj)ll rallax ). (9) 

Since the net translation of the 3-4 receiver pair is inde- 
pendent of the 12 pair, it is straightforward to arrange 
for the surviving slow oscillation to produce a net signal. 
A simple configuration (Fig. ^) puts Ax\2 = A234 = 
Ax, x 3 = -X - Ax/2, x A = —X + Ax/2. The difference 
term ff is relatively negligible and was dropped. The 
parallax terms add, producing an oscillation going like 
cos(fcAx • X(l/r% — 1/t2)). Thus with a near and a far 
source, where \/r\ — I/V2 ~ l/^i, one can measure the 
distance to the sources by interferometric parallax. Con- 
tinuing with one source after another the entire Universe 
could be mapped out with a new "cosmic distance lad- 
der." 

So far we have presented two pairs of receivers and 
signal development via "off-line" correlation calculations. 
There may be advantages to directly correlating the sig- 
nal from four receivers, cither at the amplitude or the in- 
tensity level. Let < hhhh > be the raw four-point in- 
tensity correlation. Dropping terms that oscillate rapidly, 
and with ff — * 0, a calculation gives 



< hhhh >= M + A/2 cos( 
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Here Af\ and M2 are normalizations depending on the in- 
tensity of the two sources. Note all the other terms in the 
four point correlation either vanish after statistical aver- 
aging or reduce to the two terms given in Eq. 1101 The 
remarkable cancellations in the 4-point intensity correla- 
tion show that distances can be measured in terms of a 
standard statistical description of raw data. 

Similarly, the 4-point amplitude correlation (observ- 
able at radio frequencies) consists of sums of terms with 
phases and (fjP arallax ; along with "sum-phases" of the 
form (k\ + fc 2 ) • Ax. Sum-phases in conventional long- 
baseline interferometry|j, |5| tend to cause difficulty due 
to atmospheric fluctuation effects. Intensity interferom- 
etry is much less sensitive, as positioning accuracy is set 
by the coherence time and not the wavelength . Never- 
theless both amplitude and intensity correlations contain 
^parallax _ g]g]are essentially a generalization of the 
VanCittert-Zernicke theorem^ - and both can perform 
interferometric parallax measurements. It would be pre- 
mature to assess the technological advantages of either 
scheme here. 

Earlier we mentioned that exacting measurement of 
source positions is not required. High precision of an- 
gular positions is now replaced by high precision of the 
two baselines Axi2 and AX34. How demanding are the 
tolerances, and what can be done to ameliorate them? 
An error of one part in 10 s represents a tolerable rang- 
ing error of order km on the Ax ~ 1AU baseline. 
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Proven satellite ranging techniques pfl achieve accura- 
cies superior by orders of magnitude. We note in addi- 
tion that (/>' ' might be predicted in advance by conven- 
tional means and removed by phase-shifting, heterodyn- 
ing and signal-processing strategics, cither "on-line" and 
"off-line" . Moreover, </>(°) can be observed during running 
as a useful cross-check, and being constant under transla- 
tion, can serve as a sort of absolute positioning standard. 
Another question is whether exacting timing resolution 
is needed to ensure the "same photon" enters all detec- 
tors. The answer lies in the uncertainty principle. Let 
Aw be the bandwidth of detection, under which there is a 
variation A/fi^ = Aw^ ' /uj. Suppose this error must be 
of order 10 -8 or smaller, which is the relative phase er- 
ror from position errors. Then if position and bandwidth 
errors are in tolerance, the "same wave" enters the two 
detectors in tolerably fixed phase relation. These criteria 
in fact define "same wave" and "same photon" concepts 
operationally. At optical frequencies Aw/u ~ 10~ 8 al- 
lows errors of order lOAdHz, integration over times of 
order 1 — 10/iS, which is no barrier. 

The biggest unknowns appear to be technological and 
beyond the scope of this paper. Radio astronomy is 
so sophisticated in handling noise and phase variables 
that detailed estimates must be left to experts. For 
optical intensity correlations we can nevertheless inves- 
tigate whether there might be fundamental limitations 
from photon statistics. From the uncertainty principle 
there is a relation AnA(f> > 1, where An is the fluctu- 
ation in photon number and A<p the error in measure- 
ment of a phase. Poisson statistics suggests An ~ ^Jn 
which we choose to be conservative. In effect, measure- 
ment of a phase to order one needs one photon, which we 
will estimate as ten. Meanwhile the number of photons 
detectable nfi ux decreases somewhat faster than 1/r 2 . 
When the number of photons detected falls below what is 
needed for a reliable phase we reach an upper limit on dis- 
tance measurement. If we conservatively choose sources 
separated by Ad ~ arcsec, take Ax ~ 10 4 fcm, and op- 
tical A - 5 x 1CT 4 mm, then ^ aUax „ 1, dfi ~ 10 8 , 
which is intimidating. Assuming Ax is allowed to drift, 



it must remain steady to one part in 10 8 during a sin- 
gle phase measurement. If Ace rotates once per day (10 5 
sec), each phase measurement needs 10 4 photons per sec- 
ond. For the number of photons we rescaled the cal- 
culations of Ref. 8]. We assumed a 100 m 2 aperture 
with 10% throughput efficiency, and Hubble constant 
H = 0.7 • 100km/ s/Mpc. To order of magnitude there 
are 20 photons per second in H a light for a typical bright 
galaxy with power 10 42 ergs/s at redshift z ~ 1: Insuffi- 
cient flux for a Gpc-rangc measurement. However an AU- 
scale source at 100 Mpc range is about 10 3 times more 
luminous, allowing the measurement. Moreover, approx- 
imately 10 5 such measurements can be repeated per day. 
We can see no barriers of principle. 

In summary, we have concentrated on ambitious mea- 
surements with amplitude and intensity correlations at 
the Gpc scale because they seem to us the ultimate am- 
bition. From near Earth-orbit, or perhaps with detec- 
tors fixed on Earth, it should be possible to measure 
to r > Mpc, which would be magnificent. A number 
of radio telescopes in orbit can make mutual correla- 
tions, and correlations with ground-based receivers, to 
develop stupendous resolution. The joint Japanese/US 
VSOP (VLBI Space Observatory Program) mission had 
a 21000 km orbit and an 8m telescope^!- An ambi- 
tious Russian program RADIOASTRON plans for very 
large orbits 10]. As far as we know both intensity and 
amplitude correlations can be performed with such in- 
struments. Meanwhile numerous topics in conventional 
astronomy would benefit enormously from direct mea- 
surements of distances thousands of times smaller than 
what we discuss, and correspondingly easier to imple- 
ment. While there are a host of challenging technical is- 
sues, there is every reason to believe that a new method 
of distance measurements should be feasible with current 
technology. 
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